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ATMOSPHERIC  WATER  VAPOR:  A  NEMESIS 
FOR  MILLIMETER  WAVE  PROPAGATION-2 


Hans  J.  Liebe 

U.S.  Department  of  Commerce 

National  Telecommunications  and  Information  Administration 
Institute  for  Telecommunication  Sciences^ 
Boulder,  Colorado 


Millimeter  waves  offer  an  attractive  way  of  solving  unique 
system  problems  because  of  their  ability  to  penetrate  clouds, 
smog,  dust  or  fog.  This  makes  them  a  logical  choice  over  electro- 
optical  devices  for  adverse  weather  applications.  Spectral  lines 
of  oxygen  and  water  vapor  ultimately  limit  the  atmospheric  trans¬ 
parency;  hence,  most  applications  operate  between  the  absorption 
lines  in  four  window  regions  (that  is,  24  to  48,  72  to  110,  128 
to  160,  and  200  to  260  GHz).  Observations  have  established  the 
existence  of  excess  water  vapor  absorption  (EWA) ,  which  is  most 
evident  in  these  windows.  Excess  implies  that  the  effect  is  not 
related  to  the  known  spectral  properties  of  the  water  molecule . 

EWA  is  found  to  increase  in  nontrivial  manner  with  humidity  and 
the  discrepancies  can  be  as  large  as  a  factor  of  10.  Several 
groups  (most  notably  at  the  Appleton  Laboratory,  UK),  have 
gathered  evidence  of  EWA  from  laboratory  and  field  observations 
and  brought  forward  hypotheses  to  account  for  the  data. 

Qualitative  explanations  are  based  on  the  assumption  that  water 
molecules  in  moist  air  form  larger  molecules  with  a  dimer  being 
the  first  step  in  a  series  of  stable  species.  Hydrogen  bonding, 
ion  attraction  and  attachment  of  the  polar  H20  to  foreign 
particles  (aerosol  growth)  are  the  ordering  forces  considered  in 
the  clustering  process.  An  assessment  of  the  current  EWA  picture 
will  be  given  and  avenues  of  research  attacks  are  discussed  to 
solve  the  enigma  in  the  quantitative  description  of  the  inter¬ 
action  between  millimeter  waves  and  moist  air. 

f  ^  i 
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I .  INTRODUCTION 

Atmospheric  water  in  both  vapor  and  liquid  states  is  the 
major  deterrent  to  an  unrestricted  exploitation  of  propagation 
of  millimeter  and,  more  so,  of  infrared  wavelengths.  For  most 
applications,  the  operation  of  ground-based  systems  is  limited 
to  seven  sindow  regions  W1  to  W7 ,  these  being  the  gaps  between 
molecular  absorption  lines  and  bands: 
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The  main  attraction  of  millimeter  wave  systems  is  their 
ability  to  penetrate  the  somewhat  opaque  atmosphere  (haze,  fog, 
clouds,  dust,  smoke,  light  rain)  under  circumstances  in  which 
electro-optical  and  infrared  systems  normally  fail.  Accurate  and 
detailed  knowledge  of  atmospheric  transmission  is  essential  to  an 
evaluation  of  the  advantages  of  millimeter  waves  over  the  shorter 
wavelengths. 


One  objective  of  this  paper  is  to  focus  attention  on  unsolved 
problems  in  the  construction  of  a  valid  model  for  atmospheric 
millimeter  wave  transmission,  subject  to  the  following  conditions: 
Frequency,  v  =  10  to  1000  GHz,  with  special  emphasis  on  the  EHF 
range,  V  =  30  to  300  GHz:  altitude,  h  =  0  to  30  km;  and  relative 
humidity,  RH  *  0  to  <  100%.  The  radio  engineer  uses  the  optical 
term  "clear"  for  moist,  but  unsaturated  air  (RH  <  100%) ,  and 
ignores  haze  conditions.  Although  many  constituents  contribute 
to  the  total  atmospheric  attenuation  rate,  a  (dB/km) ,  absorption 
in  the  window  is  dominated  by  water  vapor  and  is  of  greatest  con¬ 
cern  for  practical  situations.  At  present,  describing  water 
vapor  absorption  mechanisms  and  relating  them  to  measurable 
quantities  is  partly  an  empirical  matter  and  lacks  credibility  in 
a  general  sense. 

The  following  interaction  effects  between  millimeter  waves 
and  air  are  well  documented: 
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Case  References 


(a) 

Absorption  line  theory 

...  1,  4,  13,  31 

(b) 

Moist  air  refraction 

...  2,  5,  7,  47- 

•50 

(c) 

Water  vapor  absorption 

.  .  .  1,  3,  4,  7, 

14,  16, 

17, 

24,  25 

(d) 

Water  vapor  dispersion 

...  2,  4,  5,  7, 

17,  62, 

67 

(e) 

Water  vapor  continuum 

absorption  . 

...  11,  12,  18, 

21,  22, 

24, 

26,  31 

(f) 

Oxygen  absorption  .  .  . 

...  15,  17,  19, 

65 

(g)  Complete  line  parameter  com¬ 
pilation  for  all  atmo¬ 
spheric  molecules . 110 

known  as  the  AFGL  Tape, 


HITRAN  Model  and  FASCODE  .  10,  14,  25,  29 

(h)  Scattering  loss  by  particles  33,  32,  35-37 

(i)  Dielectric  loss  of  water  .  .  25,27,  32,  102-106 
Other  physical  effects,  which  have  been  taken  into  account  in 
infrared  transmission  models  (9),  (20),  and  (23)  are: 

(k)  Submicron  aerosol  particles  .  38,  41,  42 

(l)  Growth  of  submicron 


particles  in  moist  air  .  .  35-37,  34 
Of  a  more  speculative  nature  are  considerations  of 


(m)  Watery  dimer  (H20>2  spectra  .  6,  8,  21,  43,  58,  64,  66, 

69,  81,  86,  89,  109 

(n)  Water  cluster  (H_0)  ,  n>2  .  .  21,  22,  44 

/.  n 

(o)  Water  ion  activity . 28,  45 


The  water  polymer  hypotheses  have  the  inherent  weakness  that 
except  for  the  dimer  (43),  neither  definite  structural  observa¬ 
tions  have  been  made  nor  millimeter  wave  spectra  calculated.  In 
the  cases  (m)  to  <o) ,  production  and  growth  mechanisms  (for 
example,  ion-induced)  as  well  as  size  and  number  density  distri¬ 
butions  have  not  been  identified  for  the  natural  atmospheric 
environment. 
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Measured  absorption  in  the  windows  W1  to  W7  could  provide 
clues  for  the  understanding  of  which  basic  physical  mechanisms  ’ 
are  missing  in  current  models.  New  instrumentation,  such  as  the 
dispersion  spectrometer  (63,65)  and  the  saturation  hygrometer 
(107),  might  prove  helpful  in  future  investigations.  Available 
experimental  results  on  absolute  attenuation  rates,  path  trans- 
mittances,  and  thermal  emission  from  laboratory  (46-71)  and  field 
(73-101)  observations  in  cloudless  moist  air  are  generally  higher 
than  values  predicted  on  the  basis  of  molecular  absorption  alone. 
It  is  a  fact,  even  after  allowing  for  difficulties  in  measuring 
the  highly  variable  water  content  quantitatively,  that  the 
absorption  by  water  vapor  in  the  atmospheric  transmission  windows 
W1  to  W7  is  not  completely  understood. 


i 

i 


i 
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Section  XI  of  this  paper  presents  the  framework  for  a  radio 
path  modeling  scheme  of  which  details,  such  as  measurements  and 
validation,  have  been  published  (15,17,19).  The  model  serves  as 
a  basis  for  the  definition  and  discussion,  in  Section  IV,  of  the 
problem  of  excess  water  vapor  absorption  (EWA)  in  the  light  of 
experimental  data.  Section  III  is  devoted  to  a  brief  presenta¬ 
tion  of  relevant  physical  properties  of  atmospheric  water  vapor. 


II.  EHF  RADIO  PATH  MODELING 

Millimeter  waves  traveling  through  the  cloudless  (RH  £  100%) 
atmosphere  suffer  both  deterministic  and  random  variations  in 
amplitude  and  phase.  Radiation  is  absorbed  and  refracted  by 
gases  and  submicron  particles.  This  interaction  is  modeled  with 
the  intent  to  predict  the  following  frequency-dependent  propaga¬ 
tion  effects:  (a)  absorptive  loss  of  coherent  radiation, 

(b)  time  of  propagation  between  two  points,  (c)  refractive  ray¬ 
bending  and  ducting,  (d)  generation  of  incoherent  noise,  and 
(e)  scintillations  due  to  random  fluctuations  of  the  medium  in 
space  and  time.  The  array  of  mathematical  and  empirical  expres¬ 
sions  used  to  compute  these  various  effects  is  called  the 
Propagation  Model.  Such  a  parameterization  scheme  depends  fore¬ 
most  on  spectroscopic  information  about  kind  and  number  density 
of  absorber  species  and  on  their  distribution  within  the  path 
volume. 


A.  The  Propagation  Model 

Complex  refractivity  N(ppm)  is  a  convenient  macroscopic 
measure  of  the  interaction  between  millimeter  wave  radiation  and 
the  individual  absorbers  in  moist  air.  A  value  of  N  accounts  for 
the  effectiveness  and  number  density  of  a  particular  absorber 
population.  Later  these  dependences  are  formulated  in  terms  of 
measurable  quantities.  At  this  point,  N  is  assumed  to  be  known 
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and  the  basic  relations  that  determine  the  electromagnetic 
bdhavior  of  a  radio  path  are  written  down.  The  complex  refrac- 
tivity  of  a  gaseous  medium  in  ppm 

N  =  Nq  +  D (v)  +  j  N"(v)  (1) 

consists  of  three  components;  namely,  the  frequency  independent 
refractivity  N0  plus  various  spectra  of  refractive  dispersion 
D(v)  and  absorption  N"(v).  The  atmosphere  is  characterized  as  a 
linear  network  with  both  passive  and  active  properties.  The 
amplitude  and  phase  response  of  a  plane  radio  wave  traveling  the 
distance  L(km)  and  having  an  initial  field  strength  Eq  is 
described  by 


D  =  Eq  exp  (I'D  (2) 

where 

T  =  j  (2ttv/c)  (1  +  N  •  10~6)  (3) 

is  the  propagation  constant  of  the  intervening  medium,  c  being 
the  speed  of  light.  Usually,  real  and  imaginary  parts  of  N  are 
separated  and  expressed  as  the  power  attenuation  rate  in  dB/km 

a  =  0. 1820V IM (n)  (4) 

and  the  phase  delay  rate  in  radian/km 

*  =  0. 0209VRE (N)  (5) 

The  frequency  v  is  in  gigahertz  (GHz)  throughout  the  paper. 

The  attenuation  rate  a  is  the  more  familiar  quantity  in 
atmospheric  wave  propagation.  The  phase  rate  <r  must  be  consid¬ 
ered  when  it  varies  spatially  (for  example,  radar  pointing 
accuracy,  long-baseline  interferometry,  maximum  dish  size  for 
reflector  antennas,  etc.)  or  with  frequency  (for  example,  band¬ 
width  limitations  of  a  coiranunication  channel).  Characteristics 
of  short,  horizontal  radio  paths  may  be  approximated  by  an 
average  value  of  N.  More  general  path  geometries,  such  as  a 
ground-to-satellite  link,  are  treated  by  dividing  the  path  into 
segments  having  quasi-constant  N-values  and  sunning.  The  cummu- 
lative  behavior  between  the  ray  points  sj  and  S2  is  expressed  in 
dB  by  the  total  attenuation 


A 


ds 


(6) 


or  by  the  transmittance  (multiply  by  1/10  log  e  +  0.023026  to 
convert  dB  to  Np 
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x  =  exp(-0.230A)  (7)  ^ 

The  medium  becomes  transparent  when  T  approaches  one  and,  on  the 
other  hand,  opaque  for  T  -  0. 

The  total  phase  change  for  the  same  path  is  given  in  radians 
by 

<J)  =  <p  (s)  ds  (8) 


which  translates  into  the  travel-time  in  ns  of  the  wave  by 

0  =  4>/2itv  (9) 

The  path  differential  ds  is,  in  practice,  an  increment  As  over 
which  N  is  quasi-constant,  and  depends  upon  the  altitude,  h;  the 
starting  angle,  0,  from  the  zenith  in  the  case  of  a  slant  path; 
and  refractive  bending  of  the  ray  (Snell's  law)  due  to  gradients 
5*/6s  (19). 


Absorption  by  the  atmosphere  causes  emission  spectra.  Each 
unit  of  volume  maintains  thermal  equilibrium  with  its  environment 
via  collisions;  hence,  the  path  element  ds  radiates  an  equivalent 
blackbody  emission  T(s)  a(s)  ds,  which  is  reduced  by  the  trans¬ 
mittance  T (s)  along  the  path  of  observation.  The  resulting 
brightness  temperature  in  K 


tb  - 


0.230 


T (s)  a(s)  T(s)  ds 


(10) 


is  either  less  than  or  equal  to  the  ambient  temperature  T. 
Equations  (6) ,  (8) ,  and  (10)  constitute  the  key  by  which  perfor¬ 
mance  limitations  of  EHF  systems  operating  over  clear-air  propa¬ 
gation  paths  may  be  evaluated.  A  transfer  function  exhibiting 
constant  amplitue  (A  =  constant) ,  frequency-linear  phase  delay 
(D(V)  =  0] ,  and  no  noise  (TB  =  0)  implies  ideal  channel  behavior. 
A  broadband  signal  occupying  a  frequency  interval  Av  is  distorted 
by  the  deterministic  spectra  of  D(v)  and  N"(v);  in  addition,  TB 
imposes  detection  limitations.  As  a  bonus,  the  emission  spectra 
Tg(v)  afford  opportunities  to  sense  remotely  the  state  of  the 
atmosphere  by  passive  radiometric  means  (for  example,  Ref.  39 
and  40) . 


B.  The  EHF  Refractivity  of  Moist  Air 


The  physical  state  of  moist  air  is  described  by 
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dry  air  pressure  (1  kPa  =  10  mb) 
relative  inverse  temperatue  (T  in  K) 
water  vapor  partial  pressure 

Calculation  of  the  frequency-independent  refractivity  in  ppm  in 
Eq.  (1)  is  straightforward  (47-50) 

N  =  2.589  p  t  +  (41.6  t  +  2.39)  et  (12) 

o 

Water  vapor  refractivity  is  about  16  times  more  effective,  on  a 
per  molecule  basis,  than  dry  air  in  generating  propagation 
phenomena  such  as  delay,  ray  bending,  ducting,  scintillations, 
etc. 


P  ,  kPa  •> 

t  =  300/T  l  (11) 

e  ,  kPa  ' 


The  dispersion  contribution  in  ppm 
D(V)  =  I.(S  F').  -  41.6  et2  (13) 

l  l 

and  the  absorption  spectrum  in  ppm 

N"(v)  =  E  (s  F")i  +  N"  +  N"  (14) 

require  further  elaboration.  Frequency-dependent  molecular 
spectra  are  of  two  types: 

1.  Line  spectra  of  absorption  SF"  and  of  dispersive  refrac¬ 
tion  SF' ,  having  strength  S  in  units  of  kHz  and  shape  factors  F' 
and  F”  in  units  of  GHz~l;  the  sums  over  i  consider  millimeter 
wave  lines  (see  Tables  1  and  2)  of  02  (i  =  2  to  45) 2  and  H2O 

(i  =  46  to  74).  spectra  of  the  trace  gases  03,  Co,  N2O,  S02, 
NH3,  etc.,  are  neglected  (10,17,29,31,110). 

2.  Continuum  water  vapor  spectrum  due  to  far-wing 
contributions  of  very  strong  infrared  lines.  A  third  term,  NJJ, 
which  is  not  fully  understood,  was  added  to  the  absorption  to 
account  for  contributions  other  than  those  of  the  rotational 
water  vapor  line  spectrum  (see  Section  IV) . 


Common  to  each  spectroscopic  feature  is  an  intensity-against- 
frequency  distribution  function,  the  shapes  F’(V)  and  F" (v) . 

For  h  <  20  km,  the  shape  functions  are  (17) 


F*  = 


(v  -v)  +  yi 

r_° - 

L  2  2 

<vq  -v)  +  y 


+ 


(V 


+  v)  +  yi 

2  2 

+  v)  +Y 


] 


(15) 


and 


2 


The  nonresonant  oxygen  spectrum,  i  =  1  is  discussed  in 


Ref.  17. 


TABLE  I.  Data  Base  for  oi6  Spectral  Lines  in  Air  up  to  1000  GHz 
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A  NEMESIS  FOR  MILLIMETER  WAVE  PROPAGATION 


v  p  Y  -  (vQ  -  v)  i  y  -  <vo  +  v>  X  ^ 

F"  =  (  )  Q  2  2+  2  2-1 

v  (v  -  v)  +  Y  (v  +  v)  +  y 


introducing  the  spectroscopic  parameters  of  center  frequency  vQ, 
width  Y»  and  overlap  interference  I.  F'  and  F“  are  in  units  of 
1/GHz . 

For  h  >  20  km,  only  isolated  resonant  lines  are  present, 
spread  over  a  megahertz  frequency  scale.  Equations  (15)  and  (16) 
reduce  to  Lorentzian  shapes 

F'  =  (v  -  v)/((v  -  v)2  +  y23  (17) 

L  o  o 


F"  =  Y/[(v  -  V)2  +  Y2]  (18) 

L  o 

where  F£  and  F£  are  in  units  of  1/GHz.  Peak  dispersion  at 
V  =  VQ  ±  Y  and  maximum  absorption  at  V  =  VD  are  given  in  ppm  by 


N"  =  S/Y  (19) 

m 

For  h  >  40  km,  a  further  decrease  in  pressure  converts  the 
pressure-broadened  Lorentzian  into  a  Doppler-broadened  Gaussian 
shape  with  a  different  width  in  kHz  of 

Yd  =  6.20  VQ/*(mt  (20) 

where  m  is  the  molecular  weight  (e.g.,  O2,  m  =  32,  t  =  1, 

VQ  =  60  GHz:  Yd  =  65 • 8  kHz). 

The  convolution  <  f  Lorentzian  and  Gaussian  shape  functions  is 
called  the  Voigt  profile,  which  is  governed  by  the  parameter 
Y/Yd*  The  Voigt  profile  is  appropriate  when  this  ratio  falls  in 
the  range  between  =  10  and  0.1.  Numerical  evaluation  of  the 
height-dependent  complex  Voigt  function  requires  considerable 
computational  effort  (19). 


1.  Microwave  Spectrum  of  Oxygen  (O^-WS) .  The  02  molecule 
has  electronic,  vibrational,  and  rotational  energy  levels  with 
transitions  causing  spectral  lines  from  the  optical  to  the  EHF 
range.  The  EHF  lines  are  fine  structure  transitions  between 
rotational  triplet  energy  levels.  All  pertinent  information  on 
the  02-MS  is  tabulated;  that  is,  the  line  parameters  have  been 
reduced  to  five  coefficients  a^  to  aj,  which  are  listed  in 
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Table  I  together  with  the  center  frequencies  vD  (17,19). 
Dependence  on  atmospheric  properties  is  expressed  by 

• 

3 

S  =  a1  p  t  exp(a2(l  -  t)  ] 

(21) 

0.9 

Y  *  a3  (p  +  1.3e)  t: 

(22) 

I  -  a4  p  t  5 

(23) 

where  S  and  y  are  in  kHz  and  GHZ,  respectively, 
influences  the  O2-MS  through  Eq.  (22). 

Water  vapor 

2.  Microwave  Spectrum  of  water  Vapor.  From  the  AFGL  line 
parameter  compilation  (10) ,  1838  H2O  lines  were  read  up  to  31  THz 
when  a  maximum  intensity  (v  =  \>0 )  cutoff  of  2  x  10“  3  dB/km 
(300  K)  is  applied.  For  the  EHF  range,  29  lines  (see  Table  2) 
must  be  considered  explicitly  and  the  remaining  ones  are  lumped 
into  a  continuum  far-wing  contribution.  An  individual  line  is 
described  by 

Sv  =  bj  e  t3,5  exp [b2 ( 1  -  t) ) 

(24) 

Yv  =  b3(p  +  4.80  e)t0'6 

(25) 

1=0 

v 

where  Sv  and  yv  are  in  kHz  and  GHz,  respectively. 


3.  Continuum  Water  Vapor  Spectrum.  The  remaining  1809 
(that  is,  1838  -  29)  H2O  lines  make  far-wing  contributions;  they 
are  fitted  by  (17) 

N"  .  =  1.9  x  10'5  pet4  (V/30)  (26) 

v,  x 

The  continuum  absorption  N”  in  ppm  is  of  a  magnitude  and 
functional  form  that  is  similar  to  the  empirical  Gaut-Reifenstein 
expression  (11) 

N"  =  5.6  x  10'5  p  e  t3-1  (V/30)  (27) 

v 

which  was  introduced  by  Waters  (31)  and  has  proven  useful  to 
other  workers  (32,108).  The  far-wing  contributions  to  refrac- 
tivity  of  all  rotational  H2O  lines  yield  with  Eq.  (12)  in  Eq.  (13) 
the  term  41.6  et^  (3,4,5).  The  contribution  of  the  local  lines, 
i  =  46  to  74  to  this  value  is  in  ppm 

E  (SF'  (v  -*■  01)  .  =  4.75  et3 

,  1 


(28) 
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Very  accurately  measured  refractivity  can  serve  as  a  test  for  the 
low  frequency  response  of  theoretical  line  shapes  applied  to  the 
rotational  H2O  spectrum  if  one  is  sure  of  the  spectroscopic  data 
base.  A  more  detailed  discussion  of  Eqs.  (26) and  (27)  follows  in 
Section  IV. 


4.  Temperature  Behavior.  A  closer  look  at  the  spectral  line 
response  aids  in  understanding  the  temperature  dependence  of  EHF 
attenuation  and  phase  dispersion  rates  in  moist  air.  Line  inten¬ 
sities  are  sensitive  to  temperature  variations  at  various  rates. 
The  peak  line  intensities  in  ppm  (Eq.  (19))  are  described  by 


S/Y 

°2 

(a1/a3)[p/(p  +  1.  3e)  ]  t2-1exp[a2(l  -  t)  ] 

h2o 

(b1/b3)[e/(p  +  4. 8e)  ]  t2'9exp[b2<l  -  t)  ] 

The  line  intensities  are  independent  of  temperature  when 


(29) 


t  =  1.1  1.3  1.5  2 

°2 

a2  =  2.1  in  t/(t  -  1) 

2.00  1.87  1.70  1.46 

(see  Table  1) 

H2° 

b2  =  2. 9  Jin  t/(t  -  1) 

2.76  2.54  2.35  2.02 

(see  Table  2  ) 

Lines  with  values  lower  than  those  prescribed  by  Eq.  (30) 
increase  in  intensity  when  the  temperature  drops  (t  >  1)  and  vice 
versa.  A  useful  approximation  at  a  fixed  frequency  is  given  by 

a(t)  =  a (300  k)  ty  (31) 

where  the  exponent  y  is  obtained  by  fitting  model  data. 


C.  The  Radio-Path  Model 

The  path  is  assumed  to  be  in  a  spherically  stratified  atmo¬ 
sphere  in  which  each  layer  is  homogeneous.  The  integral 
expressions  (Eqs.  (6),  (8)  and  (10))  are  evaluated  by  the  layer- 
by-layer  method  (19,29)  by  using  numerical  integration  techniques 
since  closed-form  solutions  are  lacking  because  of  the  complicated 
temperature  height  profile. 

Mean  conditions  are  modeled  by  the  U.S.  Standard  Atmosphere 
1976,  and  in  situ  data  from  radiosonde  ascents  can  be  programmed 
in  directly.  Whatever  description  of  the  atmosphere  is  employed, 
model  or  data,  it  is  converted  in  the  computer  into  n-layers 
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each  having  an  assigned  set  of  values  of  p(h) ,  e (h)  £  es,  t(h). 
Values  for  a  are  calculated  by  using  Eqs.  (4),  (14),  (16)  and 
(21)  to  (25)  and  for  ip  by  applying  Eqs.  (5),  (12),  (13,)  (15)  and 
(21)  to  (25) ,  by  the  standard  line-by-line  superposition  (29)  and 
adding  the  continuum  Eq.  (27) . 

The  radio  path  is  assumed -to  be  a  ray  (that  is,  plane  wave 
case),  and  the  general  procedure  is  to  calculate  the  local  N(h) 

(1)  for  the  programmed  altitude  grid  reporting  p,  e,  t  and  to 
store  NQ,  D (v) ,  and  N" (v)  separately  on  file.  The  ray  starts 
at  a  surface  height  hQ  with  an  elevation  angle  0  and  is  guided 
through  the  inhomogeneous  medium  until  it  reaches  the  intended 
fianl  height  h^.  Refraction  is  most  pronounced  at  the  lowest 
heights  and  causes  substantial  ray  bending  when  0  is  close  to  0° 
(tangential  path);  thus,  many  fine  steps  are  required  initially. 

Modeling  is  introduced  for  the  purpose  of  predicting  the  mean 
of  EHF  propagation  effects  and  the  limits  of  their  variability 
from  readily  available  climatological  data  bases.  In  order  to 
apply  such  a  tool  in  an  optimum  manner,  it  is  important  to  con¬ 
sider  the  uncertainties  in  the  spectroscopic  data  and  the 
limitations  imposed  by  the  numerical  calculation  procedure. 


D.  Typical  Examples 

The  examples  presented  in  Figs.  1  to  7  almost  speak  for 
themselves.  Figure  1  is  a  presentation  of  the  complete  rota¬ 
tional  water  vapor  spectrum  for  a  typical  sea-level  condition. 

It  was  calculated  by  using  the  AFGL  Data  Tape  (10)  rather  than 
the  EHF  model  given  in  sections  B  and  C,  and  serves  to  demon¬ 
strate  the  very  strong  (over  6  orders  of  magnitude)  absorption 
due  to  water  vapor  in  the  far-infrared  that  allows  on  the  low 
frequency  end  just  a  few  eindow  ranges  for  transmission. 

A  more  refined  modeling  result  is  depicted  in  Fig.  2.  It 
displays  at  higher  altitudes  (h  =  16  km)  spectral  signatures  of 
the  trace  gases  O2,  CO,  N2O,  which  are  not  included  in  the 
described  model.  At  tropospheric  heights  (h  =  0  and  4  km),  only 
lines  of  H2O  and  O2  (see  markings  at  bottom  of  figure)  are 
important.  The  figure  gives,  in  essence,  a  picture  of  atmo¬ 
spheric  molecular  absorption  in  the  millimeter  and  submillimeter 
wavelength  range.  These  computations  by  Burch  and  Clough  (24) 
agree  well  with  reported  data  and  the  presented  model  (Eqs.  (27) 
and  (4))  of  the  water  vapor  continuum  spectrum  labeled  ax.  The 
attenuation  rate  spans  seven  orders  of  magnitude  over  the 
frequency  range  v  =  100  to  1000  GHz  and  altitudes  h  =  0  to  16  km 
where  >99.5%  of  the  atmospheric  water  vapor  is  contained.  The 
transparency  in  the  window  ranges  W2  to  W5,  which  are  valleys 
between  the  absorption  line  peaks,  is  dominated  by  the  water 
vapor  continuum  ax.  Figure  2  affirms  the  dominant  role  that 
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FIGURE  1.  Attenuation  rate  a  at  sea  level  (h  =  0  km,  e  = 
1.015  kPa)  over  the  frequency  range,  V  =  30  GHz  to  30  THz  display¬ 
ing  the  envelope  for  maximum/minimum  values  of  the  rotational 
water  vapor  spectrum.  Calculated  by  using  the  AFGL  Tape  (10)  . 
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FIGURE  2.  Attenuation  rate  a  at  the  altitudes ,  h  =0,  4 ,  and 
16  km  over  the  frequency  range,  v  =  100  GHz  to  1000  GHz.  The 
following  atmospheric  conditions  and  trace  molecular  number 
densities  M(m~3)  were  used  in  the  calculation  (24) : 
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FIGURE  5.  Zenith  attenuation  A  through  the  first  30  km  of  the 
U.S.  Standard  Atmosphere  1976  for  dry  (RH  =  5%)  ,  moderate  (50%), 
and  humid  (100%)  air  over  the  frequency  range,  V  =  30  GHz  to  300 
GHz  (EHF) .  The  relative  humidity  RH  was  assumed  to  be  constant 
for  h  =  0  km  to  8  km  (19).  The  experimental  data  are  from  Ref. 

80  and  Table  5. 


Frequency ,  GHz 

FIGURE  6.  Zenith  phase  dispersion  =  <t>  -  4>0  (Eq.  (18))  for 
the  same  conditions  specified  in  Fig.  5.  The  frequency- 
independent  delay  time  due  to  <t>0lN0(h)]  is  (Eq.  (9)): 
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FIGURE  7.  Zenith  brightness  temperature  ( atmospheric  noise ) 

TB  downwelling  to  ground  from  a  dry  (w  -  0)  and  a  moist  ( w  =  2  cm, 
Eg.  (33))  air  mass  (31). 

water  vapor  absorption  plays  in  atmospheric  millimeter  wave 
propagation,  even  at  modest  humidities  (RH  =  46%  for  h  “  0) . 

A  more  detailed  picture  of  molecular  attenuation  appears  in 
Fig.  3.  The  oxygen  microwave  spectrum  dominates  in  the  ranges 
50  to  70  GHz  and  115  to  123  GHz.  At  levels  close  to  the  surface, 
the  60-GHz  lines  are  merged  into  an  unstructured  band  shape,  the 
maximum  intensity  of  which  is  pressure-proportional,  until  the 
lines  separate  (h  >  15  km) .  Above  15  km,  the  shielding  effect 
breaks  down  and  radio  channels  with  up  to  400  MHz  bandwidth  can 
be  accommodated  between  the  lines.  If  there  is  a  line  close  to 
the  frequency  of  interest,  isolated  line  behavior  takes  over.  In 
the  frequency-agile  applications,  it  is  possible  to  tune  to  a 
more  or  less  constant  shielding  factor  (for  example,  0.5  dB/km) 
over  the  height  range  h  =  0  to  «  30  km.  Above  h  >  30  km,  Zeeman 
splitting  has  to  be  considered  in  the  frequency  intervals 
v0  ±  5  MHz  until  the  02  lines  vanish  above  h  >  100  km  (19) . 

It  is  encouraging  to  the  author  to  see  experimentally  verified 
results  on  millimeter  wave  water  vapor  refractivity  expressed  by 
N0  +  D (v) ,  an  example  of  which  is  presented  in  Fig.  4.  After  many 
careful  measurements  of  microwave  refractivity  in  the  fifties  and 
sixties  (2,3,5,47-50,56),  this  marks  a  new  start  in  a  heretofore 
inaccessible  frequency  range.  In  comparison  with  conventional 
absorption  spectroscopy,  refractive  dispersion  measurements  are 
superior  for  absolute  intensity  studies  (15,57,65).  An  important 
result  of  Fig.  4  is  the  fact  that  the  much  stronger  lines  beyond 
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1000  GHz  determine,  by  means  of  their  far-wing  intensity,  most 
of  the  microwave  refractivity  N0  (see  Eq.  (28)). 


E.  Cumulative  Behavior 

A  standard  example  for  an  inhomogeneous  medium  is  the  one-way 
zenith  response  through  the  U.S.  Standard  Atmosphere  1976  (19). 
Figures  5  and  6  display  the  cumulative  attenuation  A  and  cumula¬ 
tive  phase  dispersion  A <p,  for  which  48  height  levels  up  to  h  =  30 
km  are  summed.  Except  for  the  vicinity  (vQ  ±  10  MHz)  of  the 
O2-MS  lines,  these  curves  represent  the  one-way  zenith  path 
behavior.  Three  humidity  profiles  p(h)  were  used  to  model 
RH  =  5,  50,  and  100%  in  each  of  26  height  layers  between  h  =  0  and 
8  km.  The  RH  drops  rapidly  below  1%  above  h  =  8  km.  The 
frequency  range  55  to  65  GHz  is  opaque  (A  >  30  dB)  for  any  system 
attempting  to  look  through.  The  transmittance,  Eq.  (7)  is 
measurable  when  the  atmosphere  is  somewhat  transparent  (A  <  30  dB) , 
and  it  can  be  determined  either  by  the  absorption  of  a  signal 
coming  from  the  outside  (sun  emission,  satellite  beacon)  or  by  the 
thermal  emission  originiating  predominately  in  8  km  thick  layers. 
Frequencies  around  57  and  63  GHz  yield  the  maximum  and  minimum 
phase  dispersion  for  the  60-GHz  band,  which  is  almost  independent 
of  the  water  vapor  content.  Actually,  the  foreign-gas  broadening 
of  the  02-MS  by  water  vapor  (see  Eq.  (22))  reduces  the  attenuation 
slightly  with  increasing  vapor  pressure  (see  Fig.  5:  A(100%)  - 
A ( 5%) ) .  Comparison  of  the  total  attenuation  A  with  measurements 
(see  Table  5)  yields  reasonably  good  agreement  when  the  empirical 
assumptions  formulated  by  Eq.  (27)  are  used  in  the  calculation. 

The  amount  of  water  vapor  absorption  for  a  given  atmospheric 
condition  can  be  estimated,  to  some  extent,  if  the  surface  vapor 
concentration  pQ  is  known. 

The  EHF  thermal  emission  by  H20,  02,  03  and  CO  was  calculated 
by  Waters  by  using  Eqs.  (10) ,  (7)  and  (4)  and  the  absorption 
coefficients  given  in  Eq.  (31).  Calculations  are  shown  in  Fig.  7 
both  for  no  water  vapor  and  for  a  total  precipitable  water  vapor 
of  2  cm.  The  integral  is  evaluated  for  h  =  0  to  60  km.  Clear-sky 
emission  varies  primarily  with  the  amount  of  water  vapor.  The 
calculations  are  for  observations  in  the  zenith  direction  from 
ground  level.  Cosmic  radiation  of  3  K,  which  is  incident  on  the 
atmosphere  from  the  top,  has  been  added  to  Eq.  (10) . 

III.  PHYSICAL  PROPERTIES  OF  ATMOSPHERIC  WATER  VAPOR 

The  clear  atmosphere  appears  in  the  EHF  radio  path  modeling 
scheme  only  as  a  p-t-e  (Section  II,  B)  parameter  system.  The 
real  atmosphere  is  an  enormous  theater  of  diverse  but  related 
and  incessant  activities  (33).  Water,  endlessly  changing  its 
phases,  is  the  main  actor  on  this  stage.  On  a  global  scale, 
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each  year  about  1  m  of  the  ocean  depth  is  evaporated,  whereas 
the  average  water  content  of  the  atmosphere  represents  a  depth 
of  roughly  2  cm.  Water  vapor  is  rapidly  exchanged  causing  it  to 
be  a  patchy,  capricious  medium  with  parcels  (scale  sizes,  10  to 
100  m) ,  blobs  (<  0.5  m)  and  strata  (-  1  m)  of  moisture  that 
probably  account  for  a  fair  share  of  the  scatter  and  inconsis¬ 
tencies  in  millimeter-wave  propagation  data  taken  from  field 
observations.  The  total  integrated  water  vapor  and  liquid  water 
in  a  vertical  column  can  be  detected  with  high  time  resolution 
(.2.  10  s)  by  ground-based  radiometry  in  the  microwave  range  (40)  . 
In  relating  millimeter-wave  propagation  and  atmospheric  water 
vapor,  one  has  to  be  generally  content  to  make  predictions  with 
no  more  than  statistical  certainty. 

The  fact  that  the  three  phases  of  water  contribute  to  the 
weather  cycle  at  the  prevailing  conditions  (p-t-e)  is  most 
fundamentally  a  consequence  of  the  molecular  structure  of  HoO. 
Water  molecules  tend  to  associate  through  hydrogen  bonds3  having 
about  one-tenth  the  strength  of  a  molecular  bond.  Very 
schmatically,  in  the  solid  phase  four  H-bonds  form  a  rigid 
lattice;  in  the  liquid  phase,  on  the  average,  two  H-bonds 
amalgamate  chainlike  links;  and  in  the  gaseous  phase,  a  chance 
exists  to  associate  singly  H-bonded  molecules  to  dimers  (21,43, 
44)  or  aggregate  into  clusters  of  preferred  sizes  (that  is,  10  to 
50  molecules)  (22,28)  under  the  influence  of  ion-activity  (45). 

If  further  notice  is  taken  that  atmospheric  air  is  never  free 
from  invisible  particles  having  a  variety  of  origins,  chemical 
compositions,  sizes,  and  affections  for  water  vapor,  then  the 
radio  engineer  will  almost  despair  at  the  prospects  of  ever 
putting  order  in  the  atmospheric  pandemonium.  The  value  of  the 
interdisciplinary  Workshop  on  Atmospheric  Water  Vapor  for  pro¬ 
viding  guidance  with  this  task  must  be  stressed. 

A.  Absolute  Humidity 

The  amount  of  water  vapor  (that  is,  absolute  humidity  p) , 
present  in  the  atmosphere  depends  upon;  (a)  evaporation  from 
surfaces;  (b)  transport  by  motions  on  various  scales,  mainly 
through  the  troposphere;  (c)  condensation- forming  clouds  and 
fog  that  causes  precipitation  fallout.  The  radio  propagation 
engineer  measures  the  amount  of  water  vapor  by  means  of  an 
average  concentration  p,  which  varies  at  sea  level  between  the 
extremes  of  0.1  (dry,  winter,  polar)  and  60  g/m3  (wet,  summer, 
tropical).  The  height  distribution  is  approximated  from  a  known 
ground  level  PQ  by 


3The  electron-rich  end  of  O  in  the  polar  molecule  H^O  attracts 
an  electron-poor  H  end  of  a  neighbor  HjO. 
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.  p(h)  =  pQ  exp(-h/2  km)  (32) 

and  the  most  realistic  measure  for  predictions  is  the  total 
precipitable  water  vapor  in  cm  (10  g/m3/km  =  1  cm)  • 

s2 

w  =  p(s)  ds  ~  0.2  p  (33) 

Jsx  ° 


The  quantity  w  is  measured,  for  example  for  a  zenith  path,  with 
a  microwave  radiometer  (39,40) .  The  variations  of  the  dry  air 
parameter  p(h)  and  t(h)  are  described  by  standard  height  profiles 
and  attention  is  called  to  the  model  of  the  U.S.  Standard 
Atmosphere  (for  example,  (33)).  The  number  density  for  the  atmo¬ 
spheric  gases  with  constant  volume  mixing  ratio  (O2,  N2 ,  CO2 , 
noble  gases;  or  referred  to  as  dry  air)  follows  directly  from  the 
p-t  combinations.  Over  the  height  range  h  =  0  to  16  km,  the 
changes  are  for  p  from  101  to  10  kPa  (1  kPa  =  10  mb)  and  for  t 
from  0.9  (333  K)  to  1.5  (200  K) .  Ozone  (O3) ,  besides  water  vapor, 
has  a  variable  mixing  ratio  which  is  modeled  separately.  Examples 
of  03  millimeter  wave  spectra  are  shown  in  Figs.  2  and  7  as 
reported  by  others. 


Fractional  fluctuations  due  to  water  vapor  paths  and  tur¬ 
bulence  are  typically  in  the  ranges 

Sp  =  +  80%  1 


fit  =  ±  1% 


(34) 


6p  =  ±  0.1%  J 

Variability  of  humidity  for  a  midlatitude  location  with  a  medium 
of  7.4  g/m3  is  diurnally  ±  1.0,  seasonally  ±  6.3,  and  locally 
±1.2  g/m3. 


2.  Molecular  Quantities.  The  molecular  world  of  sea  level 
air  is  very  empty.  The  molecular  radius  r  »  1.8  x  10"^  pm, 
average  spacing  d,  and  average  distance  JP.  between  collisions 
are  in  the  radios  r  :  d  :  £  =  1  s  18  i  320  (33) .  Most  of  the 
time,  a  molecule  is  unperturbed  by  neighboring  molecules.  The 
number  of  dry  air  molecules  per  unit  volume  (m-3)  is 

M.  =  2.415  x  1023  pt  (35) 

a 

(e.g.,  p  =  101,  t  -  1  ;  Md  =  2.44  x  1025) . 

Water  vapor  is  an  imperfect  gas.  From  thermodynamic  measure¬ 
ments,  it  is  known  that  there  are  slightly  more  HjO  molecules 
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per  unit  volume  than  predicted  by  the  ideal  gas  law.  The  correc¬ 
tion  is  made  by  introducing  the  second  virial  coefficient  B(t) . 
The  molecular  number  density  of  H2O  for  a  given  vapor  concentra¬ 
tion  (in  g/m-*) 

p  =  7.219  e  t  (36) 

-3 

follows  from  the  relationship  in  m 

Mv  =  [(2.989/p)  lO-23  +  B (t) ] _1  *  3.346  X  1022  (37) 

23 

(e.g.,  p  =  10  :  Mv  =  3.35  x  10  ). 

Very  few  values  for  B(t)  have  been  reported  (44,69): 
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The  deviation  e  led  to  the  postulation  of  a  dimer  molecule 
(H2O) j*  The  molecular  structure  (43)  and  the  millimeter  wave 
spectrum  of  the  dimer  (6,8,69)  are  well  established:  however, 
its  number  density  under  tropospheric  conditions  is  not  known. 

The  dimer  number  density  is  expected  to  depend  strongly  on 
temperature  since  the  hydrogen  bond  strength  is  rather  weak. 

Data  which  are  suspected  to  be  caused  by  a  dimeric  effect  can  be 
analyzed  for  their  percent  change-per-degree  or  for  their  power 
law  (ty)  dependence  to  allow  a  positive  identification: 
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Bond  strength  is  expressed  in  units  of  eV/molecule  or  kcal/mole 
and  the  value  expected  for  a  hydrogen  bond  lies  between  3  and  5 
kcal/mole.  The  dimer  number  density  in  m-3  was  proposed  to 
follow  an  expression  (6) 

Mp  =  M2  •  k (t)  (39) 

An  approximate  equation  in  m3  for  k(t)  was  derived  from  expres¬ 
sions  given  by  Bohlander  (69) 

k  Z  1.8  x  10"27  t5-6  (40) 

which  when  combined  with  Eqs.  (37)  and  (39) ,  yields  in  m 

Md  s  2.0  x  1018  p2  t5*6  (41) 

and  for  the  fractional  dimer  concentration 

Mq/Mv  *  6.1  x  10  5  p  t5'6  (42) 

(for  example,  p  =  20,  t  =  1  :  MD/MV  ~  x  30  )• 

Other  estimates  on  the  fractional  dimer  concentration  Mp/Mv 
are  reported  in  Ref.  (21):  0.4  x  10“  3  for  p  =  20  assuming 
5.4  kcal/mole;  and  in  Ref.  (44):  1.9  x  10”3  at  t  =  1  and  0.95  x 
10“  3  at  t  =  1.11,  when  extrapolated  from  saturated  water  vapor 
data  taken  over  the  range  t  =  0.77  to  0.84.  These  ratios 
Mp/^v  are  dose  to  values  of  e  defined  in  Eq.  (38). 


B.  Relative  Humidity 

The  examples  of  atmospheric  millimeter  wave  propagation 
(Figs.  1  to  7)  showed  clearly  the  dominant  influence  of  water 
vapor.  Some  difficulties  in  predicting  this  influence  obviously 
are  related  to  the  fact  that  water  vapor  is  a  vapor  and  not  a 
gas.  The  first  consequence  of  the  vapor  state  is  that  a  maximum 
concentration  pg  cannot  be  exceeded.  A  balance  exists  between 
two  states  called  the  saturation  point.  At  saturation,  the  rate 
at  which  molecules  evaporate  from  a  plane  surface  (liquid  or 
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solid)  equals  the  rate  of  incoming  condensing  molecules.  The 
saturation  concentration  over  water  for  atmospheric  conditions  i's 
fitted  by 


p  <  p  =  17.39  t6  10‘10-9-8^> 
—  s 


(43) 


in  g/m3  and  is  in  a  very  rough  approximation  simply  in  g/m3 


p  Z  25  t 
s 


•17 


(44) 


Equation  (43)  is  programmed  into  the  radio  path  model  (Section  IX) 
which  uses  RH(h)  information  to  calculate  with  t(h)  the  values  e 
or  p. 

The  saturation  concentration  defines  relative  humidity 
17  18 

f  =  p/Ps  ~  0.04  t  =  0.28  et  <  l  (45) 


or,  in  the  more  familiar  percent  notation,  RH  -  lOOf  <.  100%.  The 
maximum  vapor  concentration  varies  over  a  wide  range  as  a  function 
of  temperature. 


T,  K 
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t~17 
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Pg,  Eq.  43 

Ps*  Eq.  44 
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273. 16b 
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4.850 

4.840 
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250.00 

1.200 
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0.822 

0.822 

1 

aSmithonian  Hydrometric  Table  (see  (33)). 
^Triple  Point 


The  interdependences  between  RH,  ft  w,  and  T  are  shown  in 
Fig.  8.  The  humidity  condition  RH  =  100%  is  a  delicate  balance 
point  for  phase  changes.  A  drop  in  temperature  by  1  K  decreases 
the  water  vapor  density  of  saturated  air  by  about  6%.  Energy  is 
released  during  condensation  which  was  stored  in  the  random 
motion  of  H2O  molecules.  One  cubic  meter  of  a  forming  cloud 
that  converts  1  g  of  vapor  into  the  liquid  phase  releases  2.5  kJ. 

The  latent  heat  release  can  generate  electric  fields. 
Depending  upon  the  suddenness,  the  amount  of  cooling,  and  the 
water  vapor  supply,  more  or  less  violent  updrafts  (such  as 
cumulonimbus  clouds  in  thunderstorms)  feed  a  turbulent  air 
motion.  Small  aerosols  are  lofted  into  the  upper  region  of  a 
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FIGURE  8.  The  amount  of  precipitable  water  vapor  w  in  a 
radio  path  of  length,  L  =  1  km  and  the  corresponding  homogeneous 
water  vapor  concentration  p,  both  as  a  function  of  relative 
humidity  RH  for  various  temperatures  T.  The  temperature  depen¬ 
dence  of  saturation  water  vapor  pressure  es  and  maximum  vapor 
concentration  ps  (RH  =  100%)  (33)  is: 


T ,  C 

es , kPa 

P s'Q/™3 

-30 

0.051 

0.453 

-20 

0.125 

1.070 

-10 

0.286 

2.360 

0 

0.611 

4.840 

10 

1.227 

9.390 

20 

2.337 

17.270 

30 

4.243 

30.310 

40 

7.  378 

51.020 

The  broken  lines  indicate  schematically  the  range  where  water 
uptake  by  aerosol  takes  place  (Fig.  9) . 


forming  cloud  while  the  larger,  heavier  ones  remain  suspended  at 
lower  levels.  The  small  aerosols  carry  a  positive  charge,  the 
lower  levels  take  on  a  negative  charge.  The  charge  separation 
generates  high  electric  field  strengths  and  lightning  discharges 
occur  when  a  value  of  about  30  MV/m  is  reached.  Even  before, 
additional  ions  must  be  produced,  which  are  suspected  to  have 
catalytic  influence  on  the  formation  of  homomolecular  cluster 
(H20)n  with  n  =  10  to  50  (28). 


C.  Submicron  Hydrometeors 

Atmospheric  air  is  never  free  from  invisible  particles 
(aerosols)  having  a  wide  variety  of  origins,  sizes,  and  chemical 
composition,  and  most  importantly,  having  the  ability  to  convert 
water  vapor  into  submicron  hydrometeors.  To  account  for  aerosol 
activity,  one  must  follow  the  evolution  of  particle  size  spectra 


'IGURE  9 .  Scheme 
equivalent  sphere 
drop  sizes,  r  - 
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log-normal  distribution  (23,43).  A  major  process  for  the  forma¬ 
tion  of  the  N-mode  is  the  trace  gas-to-particle  conversion. 
Sulphur  dioxide,  for  example,  nucleates  with  water  vapor  to  form 
sulfuric  acid  primary  particles  in  large  numbers  (as  high  as 
10^5  m-3) ,  which  are  rapidly  (within  ms)  coagulated  by  Brownian 
motion  into  the  N-mode  (34).  This  mode  may  contain  sizable  mass 
concentrations,  pA(N)  ~  20  pg/m3  (38),  but  more  to  the  point,  it 
presents  an  enormous  surface  area  to  the  vapor  phase  and 
enables  the  aerosol  to  respond  quickly  to  changes  in  relative 
humidity  RH.  Particle  number  density  MA  and  total  surface  area 
SA  for  1  mg  of  water  distributed  in  spherical  droplets  of  radius 
r  with  1  m3  are  related  as  follows 


r,  pm 

M. ,  molec/m3 

A 

_  2 

S^,  cm 

N 

0.01 

2.4  X  1014 

3000 

A 

0.10 

2.4  x  1011 

300 

C- 

1.00 

2.4  x  108 

30 

(46) 

The  mixture  of  gas  and  suspended  particles  is  called  an  aerosol. 

Solution  droplets,  such  as  H2SO4,  are  highly  hygroscopic; 
that  is,  they  adapt  their  size  by  taking  up  water  if  the  RH  in 
the  ambient  air  is  larger  than  the  equilibrium  RH  over  the  drop¬ 
let's  surface.  In  a  reverse  situation,  their  size  will  shrink 
because  of  evaporation.  Other  hygroscopic  agents  are  salt 
particles  found  in  maritime  (NaCl)  as  well  as  in  urban  [for 
example,  (NH4>2S04]  environments.  These  crystals  undergo  a 
sudden  phase  transition  to  become  solution  droplets  at  critical 
RH  values  (RH  =  76%  and  80%  for  NaCl  and  (NH4)2S04,  respectively). 
The  aqueous  particles,  or  so-called  hydrometeors,  collect  in 
stable,  RH-dependent  sizes  predominantly  in  the  A-mode.  A  third 
type  of  aerosol  particle  is  relatively  passive  dust  grains  in 
sizes  commensurate  with  the  A  and  C  modes.  The  solid  matter 
accumulates  a  film  of  water  on  its  surface.  All  the  particle 
population  is  able  to  aggregate  water  (CN)  but  the  small, 
hygroscopic  parts  of  the  CN  (CCN)  play  the  dominant  role. 

The  total  aerosol  mass  concentration  in  air  with  a  humidity 
of  P  =  1  g/m3  can  range  from  PA  =  1  (clean  air)  to  400  pg/m3 
(polluted).  At  RH  =  50%,  an  average  mass  loading  of  about  0.1 
mg/m3  has  been  deduced  from  data  collected  at  291  locations  in 
the  U.S.  (38).  Values  between  0.5  (clean  air)  and  3  mg/m3 
(polluted)  are  reported  for  in  cloud  free  aerosol  (42) . 

Water  vapor-to-liquid  conversion  becomes  effective  when  RH 
exceeds  80%.  An  aerosol  population  can  triple,  even  quadruple 
its  dry  state  size,  and  thus  lead  to  more  than  a  hundred-fold 
increase  in  mass  concentration.  Since  the  average  size  stays 
below  1  pm,  a  more  or  less  invisible  cloud  exists,  optically 
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categorized  in  the  atmosphere  as  haze.  Optically  active  sizes 
(>  0.5  ym)  in  sufficient  number  are  only  produced  if  the  limit 
RH  =  100%  is  exceeded.  In  that  case,  the  available  amount  of 
water  (for  example,  that  due  to  a  sudden  cooling  of  saturated  air) 
is  shared  either  by  many  (>  108)  active  particles  (CCN)  and  small 
(r  =  8  ym) ,  uniform  droplets  form  as  fog  or  cloud,  or  by  very  few 
(<  108)  and  large  (r  >  100  ym) ,  destablizing  droplets  that  pro¬ 
duce  rain. 

The  water  uptake  of  aerosol  has  a  theoretical  foundation  in 
the  relationship  between  the  ambient  RH  and  the  equilibrium 
radius  of  a  particle.  Hanel  provided  empirical  data  on  the 
hydrometric  activity  of  aerosol  (35) .  He  defined  a  mass  growth 
factor  g(f  =  RH)  with  respect  to  the  dry  mass  concentration 
and  measured  in  the  laboratory  the  equilibrium  growth  behavior 
of  typical  air  samples.  Two  examples  are  depicted  in  Fig.  10. 

The  accretion  of  water  to  particles  is  substantial  when  RH 
exceeds  80%  and  gives  them  almost  pure  water  properties;  that  is, 
PA(f  >0.8)  =  pw.  Two  approximate  expressions  for  g(f)  can  be 
formulated  (35) 


FIGURE  10.  Water  uptake  [ mass  growth  factor  g(f  =  RH)  ]  of 
two  representative  aerosol  samples  t maritime :  Had ;  Continental: 
(NHj) 2SO4  +  Si02]  as  a  function  of  increasing  and  decreasing 
relative  humidity  RH  (35). 
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g (f )  =  4f2  (f  <  0.8)  (47) 

and 

g(f)  :  l/(l-f)  (f  <  0.95)  (48) 

if  differences  and  hystereses  in  the  growth  behavior  due  to  the 
chemical  makeup  are  disregarded.  A  detailed  discussion  of  models 
for  growth  factors  can  be  found  in  Refs.  36  and  37. 

By  referring  to  radio  path  modeling,  it  is  possible,  at  this 
point,  to  formulate, a  path-integrated  liquid  water  content  in  mm 

W  =  |  ds  =  |  p°  g(f)  ds  ~  |  [po/(l-f)  ds  (49) 

The  water  content  W  increases  at  the  expense  of  the  precipitable 
water  vapor  w  defined  by  Eq.  (33);  hence,  in  cm 

w(f)  =  w(f=0)  -  W  (50) 

The  amount  of  liquid  water  in  moist  air  (RH  £  100%)  is  diffici.lt 
to  measure.  Several  indirect  methods  have  been  developed.  Three 
more  promising  techniques  are 

1.  The  mass  concentration  of  a  dried  air  sample  is 
determined  by  micro-weighing  and  available  (35)  empirical  growth 
factors  g(f)  are  applied.  For  example,  at  RH  =  99%,  the  values 
for  g  range  between  65  and  120  and  yield  W  =  0.013  to  0.024  mm 
when  P%  =  0.2  mg/m3  over  a  path  length  L  =  1  km. 

2.  A  passive,  radiometric  microwave  method  measures  at  two 
(or  more)  difficult  frequencies  the  sky  noise  Tg  (see  Eq.  (10) 
and  Refs.  39  and  40) .  The  brightness  temperature  TB  is  converted 
into  zenith  attenuation  A  (dB)  by  means  of  the  spectroscopic  data 
base.  Guiraud  et  al.  (40)  perfected  the  technique,  which  uses  an 
instrument  operating  at  20.6  and  31.65  GHz.  The  retrieval  algo¬ 
rithms  for  the  simultaneous  determination  of  water  vapor  and 
water  liquid  content  are  adjusted  to  take  into  account  a  priori 
statistics  (that  is,  a  representative  radiosonde  profile) .  For 
example,  the  climatology  of  Oklahoma  City  yielded  (40)  in  cm 

w  =  6.22  A (20)  -  2.72  A (31)  -  0.01 
and  in  mm  (51) 

W  =  1.30  A ( 31 )  -  0.53  A(20)  -  0.001 

The  frequency  20.6  GHz  is  close  to  the  22  GHz  vapor  line  (Table 
2)  and  responds  mainly  to  vapor;  the  frequency  31.65  GHz  is 
remote  from  the  line  and  responds  to  the  dielectric  loss  of  water. 
This  instrument  provides  around-the-clock  monitoring  of  w  and  W. 
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A  typical  record  over  7  consecutive  days  (Oklahoma  City,  April, 
h  =  360  m,  locking  at  zenith)  displayed  for  each  24-hour  period- 
the  following  extremes  (40) : 


Day _ w,  cm _ W,mm 


1 

2.5 

to 

1 

0.01  to 

0.02 

2 

1 

to 

2.6 

0.02  to 

0.5 

3 

2.4 

to 

0.8 

0.02  to 

2 

4 

1 

to 

2 

0.03  to 

1.5 

5 

0.5 

to 

1.2 

0.1 

6 

0.8 

to 

1.3 

0.1 

7 

2.4 

to 

4.2 

0.1  to 

1.5 

The  point  of  this  example  is  that  the  low  limits  of  W  usually 
exceed  the  values  that  are  obtained  when  using  an  average  dry 
aerosol  mass  concentration  and  assumed  g(f)-  factors.  The  high 
limits  of  W  are  most  likely  due  to  visible  clouds  since  their 
size  range  (r  =  7  to  100  pm)  contains  the  bulk  of  the  water. 
Unfortunately,  the  data  were  not  accompanied  by  notes  providing 
information  or  meteorological  or  optical  conditions.  A  similar, 
satellite-borne  radiometer  yields,  over  ocean  surfaces  with  well- 
defined  emission  temperatures,  the  same  information  on  w  and  W  on 
a  global  scale  (39). 

3.  The  measurement  of  optical  transmission  gives  insight 
into  the  average  aerosol  state.  Both  the  optical  (A  =  0.55  pm) 
attenuation  rate  dp  and  the  visibility  V  =  20  dB/d^  are  related 
to  the  liquid  water  content  pw.  In  general,  Mie  scattering 
theory  has  to  be  applied  in  order  to  predict  values  of  dp  based 
on  available  size  distributions  and  the  complex  refractive  index 
of  the  particle  material  (23,33,35). 

For  the  purpose  of  this  paper,  it  is  surficient  to  bypass  the 
elaborate,  lengthy  calculations  and  discuss  the  main  features 
with  the  help  of  a  simple  empirical  expression  (27,37) 

a  ~  300  p*"*'^  or  V  x  0.067/p^'^  (52) 

p  w  w 

where  dp  is  in  dB/km  and  V  is  in  km.  The  relationships  between 

dp,  V,  and  the  path-averaged  water  content  j5w  are  approximately: 
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,  3 

p  .  ,  g/m  . . 

< 

0.001 

0.01 

0.05 

0.1 

0.5 

,  dB/km. 

< 

5 

19 

50 

75 

198 

V  ,  km  . . . 

> 

4 

1. 1* [1. 5] 

0. 4  * [0. 5] 

0.27* [0. 3] 

0. 10* [0.08] 

Clear 

Haze  ,  RH 

<  100% 

Fog,  RH  > 

100% 

* [  ]  =  experimental  value 

Frequently,  haze  and  fog  conditions  are  described  merely  by 
stating  the  measured  value  of  the  visibility  V*  which,  however, 
depends  on  the  definition  of  the  threshold  value  (20  dB  for  L  = 

1  km  in  Eq.  (52)].  Nilsson  performed  detailed  Mie  calculations 
on  five  different  aerosol  models  with  the  result  that  visibility 
V  (16  dB  threshold)  and  relative  humidity  RH  are  related  directly 
(23) .  The  range  of  RH  for  these  models  is 


V  , 

km 

50 

15 

5 

2 

RH, 

% 

20  to  73 

84  to  95 

96  to  99 

98.7  to  99.7 

In  summary,  the  salient  points  of  this  section  are: 

1.  The  physical  foundations  for  path-averaged  water  vapor  w 
(Eq.  33)  and  liquid  water  content  W  (Eq.  49)  were  traced  in 
measurable  quantities.  Both  are  expected  to  be  the  main  source 
for  millimeter  wave  absorption  in  moist  air. 

2.  The  number  density  M  of  the  major  molecular  absorbers  was 
given  for  O2  in  Eq.  (35),  for  H2O  in  Eq.  (37),  and  for  (H20)2 

in  Eq.  (41) . 

3.  The  water-vapor-to-water-droplet  conversion  in  atmospheric 
aerosol  and  the  associated  scale  problems  were  discussed  to  aid 

in  the  comprehension  of  the  stages  that  a  water  molecule  has  to 
undergo  before  precipitating  out  of  the  atmosphere.  The  key  role 
of  relative  humidity  RH  (Eqs.  43  to  48)  in  this  process  was 
underlined. 

4.  The  interdependence  between  the  optical  properties  (Eq. 

52)  of  moist  air  and  the  liquid  water  content  is  not  unique;  it 
depends  strongly  on  the  aerosol  size  distribution. 


IV.  EXCESS  WATER  VAPOR  ABSORPTION— EWA 

Water  vapor  absorption  (Eq.  14)  in  the  millimeter  wave 
windows  stems  from  the  fact  that  the  intensity  of  the  local 
millimeter  wave  lines  Nj{  =  ESF"  (Table  2)  is  insignificant  and 
that  two  other  absorption  terms ,  and  NJJ ,  become  dominant . 
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The  far-wing  response  of  the  rotational  spectrum  of  H2O  beyond 
1000  GHz  (Fig.  1)  is  estimated  by  applying  the  approximation 
y  <  <  v  <  <  v0  to  the  line  shape  F"  (Eq.  16)  and  results  in 

fh  *  2yv/v^  (53) 

where  F^  is  in  units  of  1/GHz.  The  H2O  far-wing  continuum  of  one 
line  in  ppm  follows  with  Eqs.  (24),  (25),  (53),  and  (11)  as 

N”  =  SF"  =  C^U.50  ppt3-1  +  p2  t2,1)V  (54) 


where 


C*  -  °-184  bl  b3/Vo  (55) 

For  moist  air,  the  p2-term  due  to  self-broadening  is  always 
smaller  (<  20%)  than  the  foreign  gas-broadening  term  pp.  The 
far-wing  contribution  of  the  strongest  H2O  line  (Table  2: 

0o  =  2774  GHz,  =  5.14  x  10-10)  is  for  an  atmospheric  condition 
described  by  p  =  100  kPa,  t  =  1.023  (20°C) ,  ps  =  17.27  g/m3,  at 
the  frequency  V  =  300  GHz,  =  4.77  x  1CT4  ppm  or  Oy  =  0.026 
dB/km,  to  which  the  p2-contribution  is  11%.  The  line  center 
attenuation  rate,  in  comparison,  is  (%  =  7.40  x  105  dB/km.  Most 
certainly,  line  shape  theory  (Eq.  16)  is  overtaxed  when  it  is 
applied  to  predict  relative  intensities  in  the  far-wing  over 
seven  orders  of  magnitude. 

An  unspecified  term  NJJ  was  added  to  Eq.  (14)  to  account  for 
discrepancies  between  predictions  based  on  summing  N^-terms  (54) 
and  measurements.  The  only  certainty  in  the  conflicting  evidence 
for  NjJ  is  its  correlation  with  atmospheric  humidity.  The  abbre¬ 
viation  EWA  (excess  water  vapor  absorption)  is  used  to  describe 
N£.  For  simplicity's  sake,  different,  mostly  exponential  tempera¬ 
ture  functions  have  been  reduced  to  a  power  law  tY  (Eq.  31) .  In 
this  section,  the  exponents  of  physical  (px,  ty)  and  frequency 
(vz)  dependences  are  formulated  for  absorption  models  and  com¬ 
pared  with  those  for  experimental  data.  This  procedure  appears 
to  be  one  practical  way  of  identifying  a  specific  absorption 
mechanism,  especially  since  EWA  observations  are  not  over¬ 
whelmingly  consistent. 

Two  schools  of  thought  have  evolved  to  explain  EWA: 

1.  A  molecular  approach  searching  for  water  polymers  (HjO) 
and  their  spectra  in  the  atmosphere.  Sizes  of  n=2  (dimer)4 
(6,8,21,69,81,86,89)  and  n=3  (cluster)  (21,28)  are  possibilities. 

4 

Jona  already  in  1919  had  considered  a  water  dimer  to  explain 
anomalous  dielectric  water  vapor  results  (72). 
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2.  Liquid  water  uptake  by  submicron  aerosol  particles  under 
conditions  of  high  relative  humidity  (RH  >  85%)  (9,20,23,25,30, 
35)  . 

Each  conjecture  is  supported  by  some  as  well  as  contradicted  by 
other  bits  and  pieces  of  experimental  evidence. 

A.  Definition  of  EWA 

Window  attenuation,  both  model  and  experimental,  is  fitted 
to  expression  in  dB/km  of  the  form 

a  =  0.182VN"  =  C  pX  ty  (V/30) Z  (56) 

where  x,  y,  z  are  the  proper  exponents  of  a  particular  absorption 
model. 


The  AFGL  compilation  lists  38  350  H2O  lines  from  20  GHz  to 
331  THz  (10)  ordered  in  seven  bands  (14) ,  of  which  1809  rota¬ 
tional  lines  up  to  V  =  13  THz  contribute  to  the  millimeter  wave 
continuum  (17) .  A  fit  to  these  results  yielded  Eq.  (26) ,  which 
can  be  reformulated  into  Eq.  (56)  .  Magnitude  and  exponents  for 
the  line  continuum  are  given  in  dB/km/g/m^/kPa  by 

C  »  8.4  x  10~6 
v 

and  (57) 


z 


V 


2 


These  exponents  are  based  on  the  molecular  number  density  (Eq. 
37)  displaying  the  exponents  x  =  1  and  y  =  1.  With  Eqs.  (56) 
and  (57),  it  follows  for  the  sample  conditions  above  (V  =  300, 
p  =  100,  P  =  17,  t  =  1.023)  that  =  2.75  dB/km. 


The  dimer  spectrum  is  discussed  in  detail  in  Refs.  69  and 
109.  Absorption  should  follow  dependencies  given  by  Eq.  (57) , 
but  modified  for  the  dimer  number  denisty  (Eq.  41).  The  response 
predicted  in  this  fashion  in  dB/km/ (g/m-*)  ^  is 

Cp  =  (value  of  Table  4) 

and  (58) 


z  =  values  of  Table  4 
D 
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Aerosol  liquid  water  attenuation  is  obtained  from  published 
dielectric  data  on  bulk  water  (Fig.  13)  yielding  in  dB/km  in  the 


Rayleigh  approximation  [G ' 
32). 


Ki>^  +  (n£)S  e“  -  2n*n^]  (25,30, 


a.  =  0.82Vp  £"/[(£'  +  2)2  +  (e”) 2  J 

A  W 


(59) 


Another  method  of  calculating  the  aerosol  liquid  water  attenua¬ 
tion  assumes  that  the  medium  has  a  refractive  index  n^  -  1,  and 
that  the  attenuation  in  bulk  water  is  increased  due  to  wavelength 
shortening  (Aw  =  A0/n<i) .  This  allows  to  formulate 

a,  »  (a  /n')  (W/L)  (60) 

A  w  w 

Both  methods  have  been  applied  to  the  latest  dielectric  data  on 
water  (106)  producing  somewhat  different  values  (see  Table  4) . 
Frequency  and  temperature  dependencies,  when  approximated  using 
Eq.  (60) ,  follow  in  dB/mm  from  the  rough  data  fit 

av/nw  ~  24  <v/300>°‘65  t_6 

Now,  if  the  liquid  water  content  is  estimated  by  Eqs.  (49),  (47), 
and  (45) ,  W  a  6  x  10“3  L  p2  t34  p£,  one  obtains  for  the  expected 
aerosol  response  in  the  form  of  Eq.  (56)  that 

x  £  2 
A 

yA*  28 


z.  *  0.65 
A 


(61) 


In  summary,  millimeter  wave  window  attenuation  a  might  very 
well  be  a  combination  of  up  to  five  different  contributions: 

|  EWA  | 

a  =  a^docal  lines)  +  av(far-wings)  +  aft  +  ap  +  (cluster?) 

(62) 

A  parametric  study  of  water  vapor  concentration  (p)  and  tempera¬ 
ture  (t)  dependencies  could,  in  principle,  reveal  the  following 
behavior  in  the  attenuation  rate : 
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Multiplication  factor  with  respect 


T 

°C 

P-S’3 

g/m 

P2S. 

(g/m3)2 

to  300 

K  data 

t 

t3 

Eq.  (57) 

t8 

Eq.  (58) 

t28 

Eq.  (61) 

-20 

1.07 

1.15 

1.185 

1.664 

3.890 

116.00 

-10 

2.36 

5.57 

1.140 

1.482 

2.850 

39.00 

0 

4.84 

23.34 

1.098 

1.324 

2.110 

14.00 

10 

9.39 

88.20 

1.060 

1.191 

1.590 

5.10 

20 

17.27 

298.30 

1.023 

1.071 

1.200 

'1.90 

30 

30.31 

918.70 

0.990 

0.970 

0.923 

0.76 

Data  obeying  x  -  1  and  y  -  3  but  displaying  magnitudes  different 
from  Cv  (Eq.  57)  can  be  interpreted  as  failure  of  the  line  shape 
function  (Eq.  16)  to  predict  far-wing  intensities.  Discrepancies 
of  this  nature  will  be  most  pronounced  for  dispersion  intensities 
(Eq.  15).  Wing  data  of  D(u)  drop  more  gradually  (“:  1/v)  with 
decreasing  frequency,  and  actually  blend  into  the  well-known 
refractivity  N0  (Eq.  12)  for  V  <  100  GHz,  as  seen  in  Fig.  4. 

B.  Excess  Water  Vapor  (EWA)  Evidence  from  Laboratory  Data 

Laboratory  measurements  play  an  important  role  in  verifying 
modeling  schemes  for  EHF  properties  of  moist  air.  Generally,  all 
the  spectroscopic  parameters  (>  200)  entering  into  Eq.  (1)  should 
be  deduced  from  absolute  intensity  measurements  under  well 
controlled  v-p-t-p  conditions.  By  a  judicious  choice  of  the 
experimental  variables,  it  is  possible  to  investigate  most 
parameters  separately. 

Numerous  millimeter  wave  and  infrared  studies  of  water  vapor 
and  moist  air  have  been  reported.  References  46  to  72  are 
selected  for  their  bearing  on  the  EWA  problem.  Evidence  for  EWA 
from  these  efforts  is  summarized  in  Table  3.  It  is  not  limited 
to  absorption,  but  also  shows  up  in  dispersion  spectra  D(v)  (63) 
and  in  refractivity  N0-studies  (48,49,52,56,63).  The  experiments 
are  performed  by  various  techniques.  The  radio  path  is  simulated 
in  an  enclosure  either  for  a  single-transit  or  a  multiple- 
reflection  (resonator)  passage.  Detection  sensitivity  increases 
with  path  length.  The  transmitted  energy  can  be  a  single 
frequency,  a  frequency  pair  for  differential  measurements,  or  a 
broadband  (Fourier  transform)  signal.  Main  variables  for  an 
experiment  are  either  frequency  or  pressure,  the  latter  being 


TABLE  3.  Summary  of  Laboratory  Studies  of  Continuum 
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preferable  for  EWA  studies.  Gas  mixture  control  and  vacuum 
reference  are  laboratory  advantages  for  absolute  intensity 
studies. 

Water  vapor  is  recognized  to  be  a  medium  that  is  difficult 
to  control ,  even  in  the  laboratory ,  due  to  its  attractive  force 
toward  surfaces.  The  surface  area  of  the  laboratory  enclosure 
replaces,  in  a  way,  the  micro-surface  of  an  atmospheric  aerosol 
population.  Water  molecules  do  not  ordinarily  aggregate 
spontaneously,  but  water  vapor  becomes  liquid  water  when  wettable 
surfaces  are  present  to  retain  the  impinging  molecules.  A 
threshold  of  >  2  x  10“®  cm-^  in  the  sur face-to-volume  ratio  is 
sufficient  to  form  a  continuous  interface  between  vapor  and 
liquid  (38).  This  value  is  always  exceeded  in  atmospheric  air 
(see  Eq.  46)  as  well  as  in  a  laboratory  test  chamber.  An 
experimenter  can  select  surface  materials  that,  to  a  certain 
extent,  passivate  the  attraction  for  water  vapor.  A  systematic 
study  was  made  of  various  surface  coatings  applied  to  an  electro- 
polished  stainless  steel  (SS  304)  cavity  (3440  cm3  an(j  1265  cm2, 
S/V  =  0.37  cm“l)  evacuated  for  >  24  hours  to  10“4  torr  and 
subjected  to  pure  water  vapor,  eQ  =  2.40  kPa  at  23.0°C.  Results 
on  the  relative  amount  e/eQ  of  water-uptake  by  the  walls  and  the 
time  response  of  the  absorption  (riA<j)  and  desorption  (riDe) 
process  are  as  follows: 


Coating 

e/eQ,% 

No.b 

\d's 

nDe'S 

Teflon  FEP120  (DuPont) , 

1  coat 

O 

O 

1 

20 

180 

500 

SS  304,  electropolished 

-1.50 

30 

150 

400 

HMDS3  silanizing 

-1.35 

27 

190 

750 

Parylene  C  (Union  Carbine) 

-1.60 

32 

140 

580 

Silicone  SR240  (GE) 

-2.10 

42 

200 

550 

Teflon  REP120 ,  2  coats 

-2.90 

58 

230 

5500 

Clear  lacquer 

-4.00 

80 

300 

7200 

60  GHz  spectrometer  cell. 

-19. 5(!) 

102 

2  hrs 

5  hrs 

untreated  (63) 

(S/V  =1.28 

cm"3) 

aHexamethylidisilazane  [  (CH3)  3SH  2NH 
^Number  of  molecular  layers 

The  need  for  minimum  S/V-ratios  and  judicious  material  selection 
for  laboratory  enclosures  is  clearly  evident  upon  comparing  the 
performance  of  one  of  the  typical  spectrometer  cells  with  the 
preceding  test.  Water  vapor  surface  effects  have  been  recognized 
(52);  yet,  in  many  cases  (46,49,56),  they  were  excused  as  a 
possible  source  of  error  for  the  reported  data.  Work  close  to 
saturation  requires  a  circulating  gas-handling  system  with 
controlled  mixing  and  continuously  monitored  RH  levels  (107). 
Additional  sources  of  error  are  disturbances  in  thermodynamic 
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equilibrium:  (a)  the  vapor  heats  up  when  injected  into  an 
evacuated  cavity,  (b)  the  absorption  process  releases  heat  and  ’ 
vice  versa,  (c)  the  ambient  gas  temperature  is  lowered  during 
pump-down.  All  laboratory  results  on  spectroscopic  data 
(N0,  D,  N")  of  water  vapor  should  be  seen  in  light  of  these 
comments.  Some  selected  examples  displaying  EWA  behavior  are 
discussed  in  the  following. 


1.  Moist  Air  Studies.  Llewellyn-Jones  et  al.  (66) 
investigated  the  temperature  dependence  of  the  frequency 
V  =  213  GHz  by  studying  moist  nitrogen  (p  =  93  kPa)  over  the 
range  T  =  270  -  320  K.  The  attenuation  (Eq.  62)  for  these 
conditions  can  be  described  by  the  empirical  expression  in  dB/km 

a  =  0.22pt2  +  0.01p2t4  +  0.03p2t26  (63) 

Close  to  saturation,  the  following  values  result  from  Eq.  (63) : 


T,  K 

t 

p,  g/m3 

a, 

dB/km 

250 

1.20 

1 

0.32 

+ 

0.02 

+  3.43  =  3.8 

273 

1.10 

5 

1.33 

+ 

0.37 

+  8.94  -  10.6 

386 

1.05 

10 

2.43 

+ 

1.22 

+  10.67  =  14.3 

300 

1.00 

20 

4.40 

+ 

4.00 

+  12.00  =  20.4 

2.  Pure  Water  Vapor  Studies.  In  the  case  of  self-broadening 
an  exponent  x  =  2  (Eq.  54)  is  expected.  Mrowinski  (56)  observed 

dx  =  3  x  10-4  p2  dB/km  at  vQ  =  22.235  GHz  and  T  =  39°C,  in 

addition  to  the  line  absorption  a^.  Measurements  by  Liebe  (63) 
at  30.6  GHz  and  61.2  GHz,  T  =  300  and  325  K  yielded  in  dB/km 

a  =  4  x  10_3pt2(V/60)2  +  8  X  10*4p2t10(V/60)2'5  (64) 

The  same  experiment  gave  for  the  refractive  dispersion  (ppm) , 

D  =  N' (61. 2)  -  N' (30.6) ,  the  result 

D  =  2  x  10-3pt8  +  1.2  x  10_3p2t16  (65) 

which  is  an  average  of  the  type  of  result  exhibited  in  Fig.  11. 

The  dispersion  response  resembles  the  water  uptake  curves  (see 
Fig.  10)  published  by  Hanel  (35)  and  provided  the  impetus  for 
the  discussion  presented  in  Section  III.C. 

Bohlander  experimented  extensively  in  the  100  GHz  to  1000  GHz 
range  and  deduced  from  attenuation  data,  including  results  from 
other  investigators,  the  component  cxx/p2  shown  in  Fig.  12  and 
listed  in  Table  4  (69,  109).  He  also  calculated  (by  theoretical 
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FIGURE  11.  Dispersion  D  of  pare  water  vapor  measured  with  a 
spectrometer  cell  between  the  frequencies  61.2  and  30.6  GHz 
(63) .  A  strong  anomalous  component  is  displayed  exhibiting  a 
condensation/ evaporation  hysteresis  typical  for  water  uptake 
activity  ( see  Fig.  10). 


means)  the  envelope  of  a  rotational  dimer  spectrum,  a  fractional 
dimer  concentration  (Eq.  42)  of  10-3  being  assumed.  The  data 
fit  between  EWA  and  the  dimer  attenuation  (c/x  in  Table  4)  is 
poor;  the  spectral  shapes  do  not  match.  The  fit  is  much  improved 
if  the  aerosol  attenuation  rate  (Eqs.  59  and  49)  in  dB/km 

“a  =  (Vn4)Pw  (66) 

is  used  in  the  comparison  (a/x  in  Table  4) .  Equation  (66) 
assumes  pure  liquid  water  droplets  of  submicron  size.  Actually, 
based  on  this  fit,  it  is  possible  to  predict  a  liquid  water  con¬ 
centration  pw  that  is  needed  to  reproduce  the  EWA  data  of  Table 
4;  that  is,  in  g/m3 

Pw  a  1.1  x  10'3  p2  (67) 

At  this  point,  the  reader  might  recollect  the  arguments  that  have 
been  brought  forward  to  formulate  Eqs.  (47),  (49),  (59),  and 
(61) .  Millimeter  wave  attenuation  by  submicron  hydrometeors  is 
derived  from  Mie’s  scattering  equations  in  the  Rayleigh  approxi¬ 
mation  (33).  The  refractive  index  and  attenuation  rate  of  bulk 
water  are  presented  in  Fig.  13  (102-106) .  The  bulk  water  attenua¬ 
tion  is  higher  because  of  refractive  wavelength  shortening  in 
the  medium;  hence,  the  division  by  is  applied  for  the  gaseous 
aerosol  medium  (n'  =  1). 
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FIGURE  12.  Summary  of  reported  frequency  dependence  for 
excess  water  vapor  absorption,  compared  with  a  rotational  dimer 
band  spectrum  (Mp/Mv  =  10~ ■*)  (69)  and  the  absorption  spectrum 
o.A  of  liquid  water  in  submicron  hydrometeors  (see  Table  4). 


Frequency  v 

FIGURE  13.  Complex  refractive  index  nw  and  bulk  attenuation 
rate  aw  of  water  as  a  function  of  frequency ,  v  »  3  GHz  to 
300  THz  (composite  from  Refs.  102  to  106). 
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The  refractivity  NQ  =  40. 4p  +  6nq,  when  measured  at  40°C  and 
21.4  GHz  and  23.6  GHz  (49,56,63)  exhibited  also  the  anomalies  ' 

6Nq  : 

RH,%  <  85  91.0  94.0  97.0  98.0  99.0 

6N0,  ppm  0  0.3  0.8  1.5  1.8  2.5 

c.  EWA  Evidence  from  Field  Measurements 

Millimeter  wave  field  measurements  are  carried  out  in  three 
ways  (37-101): 

1.  attenuation  rate  a  (dB/km)  for  horizontal,  line-of-sight 
paths ; 

2.  total  zenith  attenuation  A(dB)  ;  and 

3.  sky  noise  TB(K) ,  often  employed  to  infer  the  A  value. 

Most  of  the  measurements  are  performed  at  single  frequencies; 
a  few  were  carried  out  in  a  broadband  mode  (Fourier  transform 
technique)  (86,90,93,98,101).  Several  difficulties  plague  field 
observations : 

1.  absolute  calibration; 

2.  large  scatter  (+  10%  to  ±  30%)  in  data  due  to  unspecified 
weather  along  the  path; 

3.  scarceness  of  data  at  high  humidities  (RH  >  90%); 

4 .  lack  of  simultaneous  recordings  of  the  integrated  water 
vapor  w  (Eq.  33) ; 

5.  absence  of  data  on  integrated  liquid  water  W  (Eq.  49) 
and  visibility  V  (Eq.  52)  in  cloud-free  air;  and 

6.  difficulty  of  fitting  data  empirically  to  surface-based 
meteorological  variables. 

Horizontal  and  zenith  path  data  are  summarized  in  Table  5  and 
compared  with  model  calculations  by  using  the  U.S.  Standard 
Atmosphere  (19).  As  expected,  the  window  (W1  to  W5)  absorption 
increases  with  the  water  vapor  concentration  p  or  path-integrated 
water  vapor  w.  It  was  observed  earlier  that  measurements  in  the 
100  GHz  to  117  GHz  range  revealed  considerably  higher  values  of 
water  vapor  absorption  than  were  predicted  by  H20  line  shape 
theory  (73) . 

Generally,  the  absorption  A  is  divided  into  a  dry  contribu¬ 
tion  (Ad)  and  a  wet  term  (Av) .  The  dry  term  is  caused  by  oxygen 
absorption  originating  from  a  well-known  spectroscopic  data  base 
(Section  II);  hence,  in  W1  to  W3,  it  can  serve  as  a  check  value. 
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The  wet  term  determines  the  transparency  in  all  window  ranges. 
Transparency  is  tightly  coupled  to  humidity  [that  is,  p,  w,  and 
W(RH) ]  causing  various  propagation  limitations  (for  example, 
usable  range  L  or  minimum  elevation  angle  6)  with  increasing 
amounts  of  water  vapor  in  the  path  volume.  A  typical  case  exists 
when  weather  conditions  change  from  clear  and  dry  to  cloudy  and 
wet.  The  water  vapor  absorption  problem  fades  rapidly  away  above 
h  >  3  km.  It  has  become  a  standard  practice  to  separate  the  wet 
term  into  two  components:  one  proportional  with  p  (monomer),  the 
other  with  a  (dimer)  dependence.  The  first  field  results  to 
support  the  dimer  hypothesis  were  taken  at  220  GHz  (76) . 

Comparing  these  data  with  other  available  data  in  Table  5,  one 
notices  that  in  about  half  the  cases  a  squared  (x  =  2)  water 
vapor  dependence  improves  the  fit. 

A  horizontal  path  that  operated  at  182.9  GHz,  close  to  the 
183  GHz  H2O  line  (Table  2) ,  exhibited  differences  in  the  water 
vapor  proportionality  depending  upon  clear  or  cloudy  sky  condi¬ 
tions.  The  difference,  3.9P  compared  with  4.2p  dB/km,  was 
reconciled  when  a  temperature  dependence  of  t10  was  assumed, 
which  is  about  four  times  the  t  dependence  for  the  183  GHz  line 
(Eq.  29:  t2-3)  (77).  The  same  experiment  produced  at  171  GHz  a 
relationship  between  attenuation  rate  p  and  concentration  that 
becomes  increasingly  nonlinear  when  p  >  10  g/m3.  Measurements  in 
W5  (330  GHz  to  360  GHz)  required  a  linear  (x  =  1,  y  =  2.3)  and 
squared  (x  =  2,  y  =  11)  water  vapor  term  to  fit  the  data  (89). 
Condensation  effects  were  evident  in  recent  W4  and  W5  data  (101) . 
For  the  first  time,  the  importance  of  relative  humidity  RH  was 
recognized.  An  EWA  contribution  with  x  =  2  and  y  =  16  to  30 
can  be  isolated  and  an  increase  of  ax  with  beginning  fog  condi¬ 
tions  was  measured.  Predicted  attenuation  rates  in  fog  are  in 
the  range  3  to  15  dB/km  for  V  =  100  to  300  GHz,  pw  =  1  g/m  and 
T  =  4“C  (25) .  An  independent  measurement  of  water  vapor  w  and 
liquid  W  contributions  would  be  desirable. 

Cumulative  data  taken  through  the  total  air  mantle  are  more 
difficult  to  interpret  for  their  phenomenological  origins.  In 
most  cases,  zenith  attenuation  A  is  correlated  with  the  surface 
water  vapor  concentration  P0 .  Simultaneous  measurements  of  the 
integrated  vapor  w,  either  by  microwave  (see  Eq.  51)  or  infrared 
techniques  (93)  are  rare.  The  maximum  attenuation  detectable  in 
the  EHF  range  is  A  £  30  dB.  The  dry  term  A^  of  a  slanted  radio 
path  follows  the  secant  law  Ad(9O“)/sin0,  even  for  low  angles 
(0  S  10°),  the  wet  term  Ay  increases  very  rapidly  below  0  <  10“ 
(91) .  A  tangential  path  (0  =  0)  traverses  about  38  times  the  air 
mass  at  zenith  but  perhaps  100  times  the  water  vapor  content. 

It  seems  certain  that  Ay  is  a  continuum  spectrum  (79,93,98, 
108)  and  that  earlier  observations  of  spectral  dimer  features 
(81,86,90,64)  were  instrumental  effects.  To  quote  from  Ref.  93: 
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" [EWA]  shows  no  spectral  features  and  scales  with  w. "  The 
frequency  dependence  of  the  empirical  continuum  spectrum  * 

(Eqs.  28  and  4)  was  z  -  2.  A  fit  of  the  A v  term  between  15  GHz 
and  230  GHz  (Table  5)  yields  in  dB  approximately 

Ay  -  0.013  po(v/30)1*7  (68) 

and  the  frequency  dependence  drops  to  z  =  1.22  for  window  data 
taken  between  140  GHz  and  300  GHz  (93) . 

Field  data  mostly  serve  the  practical  purpose  of  establishing 
an  operational  data  base.  Agreement  with  model  calculations 
(Section  II)  is  fair,  when  it  is  considered  that  the  model  makes 
use  of  Eq.  (27)  to  describe  the  water  vapor  continuum.  Either 
integrated  vapor  w  or  surface  concentration  p  are  useful  predic¬ 
tors  of  millimeter  wave  window  absorption  (Table  5). 

D.  EWA  Discussion 

A  water  vapor  continuum  spectrum  dominates  the  transmission 
behavior  of  atmospheric  window  ranges  (W1  to  W5:  Figs.  1  to  3). 
Three  absorption  mechanisms  have  been  parameterized  to  aid  in 
identifying  the  respective  relative  contributions  to  the  con¬ 
tinuum.  In  the  course  of  this  paper,  the  following  picture 
evolved  s 
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The  empirical  expressions  (Eqs.  62  to  68)  merely  organize 
experimental  data;  they  are  not  found  to  be  entirely  satisfactory 
to  uniquely  support  even  one  particular  absorption  scheme.  At 
present  it  is  difficult  to  relate  EWA  contributions  quantitatively 
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to  observed  variables.  Results  vary  greatly  from  one  experiment 
tfo  the  next,  even  under  controlled  laboratory  conditions.  EWA 
is  most  prevalent  when  RH  >  90%;  however,  the  best  type  of 
measurement  has  not  yet  been  undertaken. 

There  are  objections  to  both  interpretations  of  EWA.  Unclear 
problems  with  respect  to  the  dimer  are:  (a)  no  direct  physical 
evidence  in  the  atmosphere;  (b)  no  match  to  the  frequency- 
envelope  of  a  proposed  spectrum  (Table  4) ;  and  (c)  no  unique 
match  to  the  temperature  dependence.  A  promising  EWA  hypothesis 
is,  at  present,  the  vapor-to-liquid  conversion  of  water  by 
aerosol  particles  under  conditions  where  RH  ^  100%.  Here  the 
objections  lie  in  the  fact  that  invisible  or  haze  clouds  (see 
X-distribution  in  Fig.  9)  must  exist  with  liquid  water  concentra¬ 
tions  close  to  those  measured  in  fogs  or  clouds  (for  example, 

RH  =  96%,  T  =  10°C,  p  =  9,  and  pw  =  0.1  g/m3,  by  using  Eq.  (67)). 
Delogne  raises  the  point  that  for  the  case  of  an  inhomogeneous 
distribution  of  liquid  aerosol  water,  the  visibility  V  is 
improved  due  to  the  p% ’ 6-dependence  (Eq.  52)  whereas  millimeter 
waves  respond  just  to  the  total  water  W  (Eq.  49;  Ref.  30). 

Which  explanation  for  EWA  has  the  most  merit  can  only  be 
answered  by  additional,  completely  controlled  experiments  per¬ 
formed  in  the  laboratory. 


V.  CONCLUDING  REMARKS 

The  role  that  atmospheric  water  vapor  plays  in  millimeter 
wave  propagation  was  traced  by  means  of  data  from  modeling, 
laboratory,  and  field  studies.  The  assumptions  and  approxima¬ 
tions  made,  the  reasoning  used,  the  relative  importance  of 
various  parameters,  the  limitations  of  available  data,  and  the 
separation  of  water  vapor  and  liquid  water  effects  have  been 
addressed.  Molecular  absorption  due  to  oxygen  is  prominent  in 
the  45  GHz  to  125  GHz  range  (Fig.  3);  water  vapor  dominates  at 
higher  frequencies,  actually  up  to  30  THz  (Figs.  1  and  2). 
Considering  only  molecular  absorbers  (H2O  and  O2)  for  a  cloudless 
atmosphere  leads  invariably  to  discrepancies  between  predicted 
and  measured  attenuation  rates.  The  problem  is  most  apparent  in 
the  EHF  window  ranges.  A  horizontal  path  at  sea  level  can 
experience,  at  220  GHz  for  example,  the  following  varying 
attenuation  rates  (32): 


Clear  air 
Fog 

Rain  (<  10  mm/ hr) 


1.6  to  11.2  dB/km 
0.4  to  4.7  dB/km 
1  to  7  dB/km 


These  values  demonstrate  the  relative  importance  of  understanding 
the  clear  air  problem  first.  Clear  air  attenuation  is  related 
to  atmospheric  water  vapor  content;  unfortunately,  more  than 


two-thirds  of  this  contribution  is  described  only  by  empirical 
formulas  that  lack  both  physical  insight  and  general  appli-  • 
cability. 

A  reasonably  concise  model  for  the  molecular  absorption  was 
presented  in  Section  II  and  then  applied  to  identify  the  magni¬ 
tude  and  frequency  and  physical  dependencies  of  the  non-HjO 
contribution  to  the  available  data  body  reported  for  laboratory 
and  field  experiments.  The  contribution  is  called  EWA  (excess 
water  vapor  absorption)  and  it  was  found  the  EWA  exhibits  in  the 
millimeter  wave  range  a  continuum  spectrum  similar  in  frequency 
response  to  the  dielectric  loss  spectrum  of  water.  The  search 
for  liquid  water  in  clear  air  led  to  the  water-uptake  phenomenon 
of  aerosol  particles,  which  grow  rapidly  in  moist  air  with 
increasing  RH.  This  problem  has  been  addressed  in  the  infrared 
(9,20,23)  and  deserves  further  study  with  respect  to  influences 
pn  millimeter  waves. 

Most  of  the  evidence  found  in  parameterizing  EWA  data 
supports  a  condensation  phenomenon  driven  by  relative  humidity, 
RH:  (a)  the  nonlinear  pressure  dependence  (aerosol  growth 
function),  (b)  the  strong  negative  temperature  dependence,  (c) 
hysteresis  effects  when  cycling  the  RH,  (d)  the  governing 
variable  is  relative  (RH)  and  not  absolute  (p)  humidity,  and 
(e)  the  failure  of  molecular  approaches  hypothesizing  weakly 
bonded  dimers  but  neither  matching  the  spectral  response  nor  the 
temperature  dependence  and  number  concentration  necessary. 

From  a  practical  point  of  view,  one  is  looking  for  invisible 
liquid  water  content  in  the  range  W  -  0.02  to  0.2  mm  when 
RH  100%,  which  are  normally  known  to  exist  under  haze,  cloud 
and  fog  conditions.  But  clouds  and  fog  imply  only  larger, 
optically  active  (r  >  1  pm)  particle  sizes.  Should  the  condensa¬ 
tion  hypothesis  have  merit,  then  EWA  can  serve  as  a  tracer  to 
the  otherwise  not  directly  accessible  world  of  atmospheric  sub¬ 
micron  particles.  A  systematic  EWA  study  under  controlled 
laboratory  conditions  of  high  humidity  (RH  **  90%  to  99.9%)  could 
reveal  growth  functions  for  particular  aerosol  ensembles  and 
elucidate  the  various  stages  and  time  scales  that  H20  molecules 
undergo  in  moist  air  until  they  become  barely  visible  as  clusters 
of  >  10^  molecules. 


SYMBOLS  AND  ACRONYMS 
Acronyms 

AFGL  Air  Force  Geophysical  Laboratory 

CM  condensation  nuclei 

CCN  cloud  condensation  nuclei 
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D1-D4  O2-MS  doublets  (Table  1) 

EWA  excess  water  vapor  absorption 

E18  computer  short-form  for  1018 

EHF  extreme  high  frequency  range  (30  to  300  GHz) 

HITRAN  high  resolution  transmission  model  developed  by  AFGL 

H-bond  molecular  hydrogen  bond 

ID  quantum  number  identification 

IM  imaginary  part  of 

O2-MS  oxygen  microwave  spectrum 

RE  real  part  of 

W1-W7  atmospheric  transmission  windows 


Wy  mbols 
a 

ai-a5 

A 

b 

bi-b3 
B  (t) 
c 
C 
d 

D(V) 

e 

E 

f 

F' 

pn 

g(f) 

h 

i 

1 

j 

k(t) 

l 

L 

m 

M 

n 

n' 

n" 

N 

No 

N’(V) 

N"  (V) 

P 

6p 

Q 


coefficient  for  fitting  p,  Tabel  5 

O2-MS  line  coefficients,  Eqs.  (21)  to  (23)  and  Table  1 

total  attenuation,  Eq.  (6),  dB 

coefficient  for  fitting  data  <*  p2 ,  Table  5 

H2O  line  coefficients,  Eqs.  (24)  and  (25),  Table  2 

second  virial  coefficient,  Eq.  (37)  m3 

speed  of  light,  Eq.  (3) ,  km/s 

continuum  magnitude,  Eq.  (53) 

molecular  spacing.  Section  III,  m 

dispersion  spectrum,  Eq.  (13),  ppm 

water  vapor  pressure,  Eq.  (11)  ,  kpa 

field  strength,  Eq.  (2) ,  V/m 

relative  humidity  RH,  Eq.  (45),  1 

line  shape  for  D(V) ,  Eq.  (15)  ,  GHz-1 

line  shape  for  N" (V) ,  Eq.  (16) ,  GHz-1 

aerosol  mass  growth  factor,  Eq.  (47),  1 

altitude,  Eq.  (32),  km 

running  spectral  line  index,  Tables  1  and  2 
interference  coefficient,  Eq,  (23) ,  1 
=  >/-I ,  Eqs.  (1)  and  (3) 
dimer  equilibrium  factor,  Eq.  (39),  m-3 
mean  free  path  length,  Section  III,  pm 
radio  path  length,  Eq.  (2),  km 
molecular  mass,  Eq.  (20) 

molecular  number  density,  Eq.  (36),  molec/m3 

refractive  index,  Fig.  13,  1 

real  part  of  n,  Eq.  (58) ,  1 

imaginary  part  of  n.  Fig.  13,  1 

complex  refractivity,  Eq.  (1),  ppm 

frequency-independent  part  of  N,  Eq.  (12) ,  ppm 

refraction  spectrum,  Eq.  (27),  ppm 

absorption  spectrum,  Eq.  (14),  ppm 

dry  air  pressure,  Eq.  (12),  kPa 

fluctuations  of  p,  Eq.  (34),  kPa 

resonator  quality  factor,  Table  3,  1 
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aerosol  particle  radius,  Eq.  (35),  pm 
relative  humidity,  Eq.  (45) ,  % 
ray  path  length,  Eq.  (6) ,  km 
ray  path  coordinate,  Eqs.  (6)  and  (8) 
path  length  increment,  Eq.  (6) ,  km 
line  strength,  Eqs.  (21)  and  (24)  ,  kHz 
aerosol  surface  area,  Eq.  (46),  cm^ 

300/T,  normalized  inverse  T,  Eq.  (11),  1 
temperature,  Eq.  (11),  K 
brightness  temperature,  Eq.  (10),  K 
visibility,  Eq.  (52) ,  km 

air  volume  containing  water  vapor,  Section  IV,  cm3 

total  precipi table  water  vapor,  Eq.  (33)  ,  cm 

total  precipitable  liquid  water,  Eq.  (49),  mm 

concentration  exponent  (p) ,  Eq.  (53) 

temperature  exponent  (t) ,  Eq.  (53) 

frequency  exponent  (V),  Eq.  (53) 

attenuation  rate,  Eq.  (4),  dB/km 

line  width,  Eq.  (22)  ,  GHz 

Doppler  line  width,  Eq.  (20) ,  kHz 

propagation  constant,  Eq.  (3„ ,  1 

deviation  from  ideal  gas  behavior,  Eq.  (38),  1 

propagation  delay  time,  Eq.  (9) ,  ns 

time  constant  for  adsorption.  Section  IV. fl.,  s 

time  constant  for  desorption,  Section  IV, B. ,  s 

slant  path  angle.  Section  II,  degree 

optical  wavelength,  Section  III,  m 

frequency,  Eq.  (3) ,  GHz 

line  resonance  frequency,  Eq.  (15),  GHz 

water  vapor  concentration,  Eq.  (36) ,  g/m3 

path-averaged  value  of  p,  Eq.  (52) ,  g/m3 

fluctuations  of  p,  Eq.  (34),  g/m3 

transmittance ,  Eq.  (7),  1 

phase  rate,  Eq.  (5),  rad/km 

total  phase  delay,  Eq.  (8) ,  rad 

sum  of,  Eqs.  (13)  and  (14) 


Subscripts 


aerosol 

dimer 

dry 

absorption  line 

Lorentzian 

maximum 

integer,  cluster  size 
initial  value 
saturation 
water  vapor 
liquid  water 
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X  EWA 

y*  optical  wavelength  (X  =  0.55  ym) 
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